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Decomposition of the Schrtidinger Equation for 
Two Identical Particles and a Third Particle of 
Finite Mass 


A. K. Bhatia md A. Temkin 
Theoretical Division 
Goddard Space Flight Center 
National Aeronautics and Space Administration 
Greenbelt, Maryland 

An angular momentum decomposition of the Schr8dinger equa- 
tion is extended to the case noted in the title. (The case of 
three unequal mass particles is treated in an appendix. ) The 
decomposition is effected with the use of a symmetric choice of 
Euler angles, and the radial equations are given in two useful 
forms. The radial equations are shown to yield the Born-Oppen- 
heimer equations for H 2 + in the limit that the two identical 
particles approach infinite mass. Other aspects of this limit 
are discussed, and general rules which relate the total angular 
momentum states of the three -body system with the molecular 
states of H 2 + are examined* 
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In a previous review we have presented a detailed investiga- 
tion of the decomposition of the Schrttdinger equation for two iden- 
tical particles, obeying the exclusion principle, in the central 
field of an infinitely heavy center of force. This, of course, 
is the prototype of the two -electron atom. The major correction to 
this idealization is the recoil effects of the nucleus. It is 
the primary purpose of this paper ti give the extension of the 
above decomposition to include the finite mass of the previously 
assumed fixed center. (The changes necessary in the radial equa- 
tions for three unequal mass particles are given in Appendix B. ) 

The resulting decomposition is still sufficiently general to 
allow an arbitrary mass to the identical particles as well as the 
nonidentical third particle. As a rather different limit of this 
system from the previous one, one can consider the limit in which 
the two identical particles become infinitely heavy and the third 
particle assumes the electron mass. This corresponds to the H 2 
limit, and it is of considerable interest to see how the radial 
equations and the three -body symmetries go over into the equations 
and symmetries of the Born-Oppenheimer approximation in that case. 

This is particularly relevant at the present time since the 
absorption of |u~ mesons in hydrogen^ via the molecule (p - jv - p) 

can presumably only be explained by the deviations of this j*. mesic 

+ . 3 

counterpart of H 2 from the Born-Oppenheimer approximation . The 

H 2 + limit will be the subject of Section III. 


- 3 - 


II KINETIC ENERGY AND RADIAL EQUATIONS 

Let f t and f x be the coordinates of two identical particles, 
each of mass m. In an arbitrary fixed coordinate system. Let 
f * 3 be the coordinates of a third particle of mass M. The kinetic 
energy can be written: 


T 





( 2 . 1 ) 


With the introduction of the difference vectors 
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(2.2a) 
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and the center of mass coordinate R: 
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the kinetic energy becomes 


T - - i- [ V 


) - -h 7 * - nr 


>+M) 


V, 


(2.3a) 


vhere 
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■+ M 


(2.4) 


The last term in (2.3a) is clearly the kinetic energy of the 

center of mass. In any closed system the potential will be independent 

of R, so that the effect of this term will be to subtract a center 

of mass energy, E,_.> from the total energy, E , to give an effective 
UM J. 

energy E for a SchrBdinger equation with the last term in (2.3a) 3Lb$€Tit> 
Thus we can replace the kinetic energy by an effective kinetic 
energy governing the internal motion 


- 4 ( ) 
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(2.3b) 
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to be used with an effective energy 

E — E T _ (2.5) 

now differs in form from the kinetic energy of the two- 
electron fixed nucleus problem^ - by the addition of the final cross 
term. This additional term is the well-known mass polarization term 
which in most helium applications is treated in perturbation 
theory (which is quite adequate for the present experimental accuracy)^ 

We shall include the effect of this term exactly. 

It will be recalled from I that the major task in effecting 
the decomposition is to find the kinetic energy in terms of the 
Euler angles and residual coordinates in place of the particle 
coordinates. We introduce formally the 

same Euler angles as in I; namely 9 the angle between the space-fixed 
£ direction and r^ x r? a', § the angle between the space-fixed 

x direction, and z x z x|', and the angle between x' and (ra-rj). 

Now, since the form of V + V is the same as in I, the 

r x *2 

transformation of those terms in terms of the three Euler angles and 
three residual coordinates can simply be taken over from I. We 
need only consider therefore the remaining cross term, and since 
we handle it in a very similar manner, we shall be mercifully brief . 

In spherical coordinates the cross term is 


v *. • a cos 9, a L_ 
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where 
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(2.7) 
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4 = cos 1 4 k - gfl ) 

sw»>, Sin ^ 


(2.10) 


c.osi>, si'n(f ( ( 2 . 11 ) 

sin 


When the residual coordinates are taken as r x , r 2 ,© 12 . 


then the 
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transformation only involves first and second partial derivatives of 
i>, j with respect to 9, f , £ , ©! 2 . (The 

latter are referred to as ^ , •( = 1 ...4 below). The cross term 
can therefore he written 




= cose lx I. T 

at.-av^ L 

J- ( U 

^ x- 


•v> \ 


+ * ±.^{ 

4 

I I a 
^ * 1 

*** 


+ JL 

H 

)V 

+ C c * 

•8i>, ^ 


+ 

4 

I 
* - | 

I i ^ 

L **,**,_ 


4 4. 

(* 




a<f x 

4 

II l i l 

V** 


II 

4L 

* 

3*, 

»^ t 


o ? ) 3 -h, 

(2.12) 

+ o If + M L?* 





L't *■)}] 


■*- 


-V- 


12* \2t) 



- 8 - 


The results in the form of the coefficients of the derivatives 
with respect to \ expressed as functions of 6,$ , * , © 12 are 
given in Table I. Finally we can write the cross term as an operator 
whose effect on the vector spherical harmonics have been fully- 
explored in I. 
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/\ 1 and /\ 2 Rrp given in (164) and is the well-known total 
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angular momentum squared operator, (l 4 l), in terms of the Euler angles. 
The wave function is now expanded in the form (155) 


■K.+ 


(■*»»<*) + 


¥ if* i ) = 2 L t s 1 »*■*.. ®|») ^ C6j i 
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where 2L 8X6 the exchange vector spherical harmonics (154), 

the simultaneous eigenfunctions of M 2 , M and ^12- The operator 

nxr 2 

Y*. ’ Y*a. as expressed (2.15) is now such that its operation on 
these ^ 1 functions can readily be found from (165,66). Thus 

the radial equations can simply be derived as an addition of terms 
to what was derived in I. Below we shall give the equations in two 
forms; the first in terms of simultaneous coupled equations for the 
"radial" functions : 
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The numbers and B* are defined in (168), and S -wave 

operator L Q in (171). The fact that the equations are in this 
w 12 

first form are coupled for a given «. implies that we are dealing 
with symmetric and antisymmetric functions 


X -f 



€ + X * *** 


(2.16a) 


X* -t + X -ft w _ 

t C *•. > , ®IX ) = - ^ ( H > *1 ; ®iv ) 


( 2 . 16 b) 


where the upper sign refers to the totally space symmetric wave 
function (singlet) and the lower > the space antisymmetric solution 
(triplet). The (anti) symmetry means that we can confine the solution 
to half the r x - r 2 plane with an appropriate vanishing of the func- 
tion or its normal derivative along r x = r 2 (cf. Eqs. (172) and (173 ) ) - 
The second form of these equations that we shall give involves 
the "radial” coordinates r x , r 2 , r 12 and involves the asymmetric 

*t *+ „ 

functions and F e 
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where, on the rhs, © 12 is understood to he expressed as a function 
of rj_, r 2 , r 12 through the law of cosines. By virtue of (2.16) 
we have 


(A* A, > V ) = ± c-0 


l + *i 


7 C *. 


This relation enables us to write the equations in an uncoupled form 
for a given n (the various x's are of course coupled to each 


other) but in a domain covering the whole r - r plane. It is also 

I 

convenient in writing this second form to combine some of the terms 


multiplied by 2 - f-VM in (2.15) with L (given in (l8l)) to form a 

r 12 

new S-wave operator L 1 : 

r 12 






A- £ 


-I 2 )^ 


S. 


t [ i. i l * lx 

" L •‘.a. ^ 


± V 


+ ’ l >v ~*-i 


A., ♦*.. 




*,x i 


|i- 


(± 

M 


Z- "3 




We also use the quantity 
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The equations are: 


iu ia + *_£ (E - v; 3 V 

■t5 x 


- * tj > l {wu *° ■ n “ ) q? 1 * k i f/ 


- ) ‘'■i’ll. F.* - 


+ (A? - + ” Nv ) ^ 


VI +1 


ix (. + **■». - ^ \ * C I “ fo~ 1 £ ^ 

2. a. 


*,«. *T~ 2 


( -0 - X ) r n.{ \ 

^ <l v 1 L *• 1 


Ns. 


-I V 


n »+a 


Vix T»V 


€(?■*■'? h. ( *i.x ^ivi ^ -V- * t>l *x . J — ^ 

' r x r 


*S« W' +?”'-)} 


+ 64^ ^ 1 ^ X - ^ I — SoTt — £** ) ^ ^ 




rD 


2^. T (% x +*»1 - ‘‘•ix ■> ^ m + o - 
n L L ^x 


- j f / 


r 


"■ p ( 

j_ 2_ 

\ 

■> 

?■ K 
■) F / 

^.nx 


\ 


CV*' 

~ W.+ X 

- u 

)} 

+ 

^T. | \_a. £ | — Sox. ” ^ 


+■ S«V 


- O 


(2.21) 



Ill THE H 2 LIMIT 


Consider the limit in which the mass of the two identical par- 
ticles becomes infinite and the third particle retains its finite 
mass. We shall call this the H 2 + limit even though in a real H 2 + 
molecular ion the two nuclei are not infinitely heavy and the system 
should he described by the complete equations we have given in the 
previous section. Mathematically this limit is defined by 


m = M 
P 


M = m 
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lim |u = m 

M p© e 
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The kinetic energy becomes 


Lt. 
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Letting 


r = 


I (fa + £2) 


(3.3) 


fas = r x - r 2 . 


(3A) 


From (3*5) £ is obviously the vector from the midpoint of the tvo 
nuclei to the electron; ve find for the effective kinetic energy: 


t « „ . 

r 

2 w e 


(3.5) 


so that the SchrBdinger equation becomes 
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C 

t x 


e -i/;J 
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(3.6) 


This is, of course, the Born-Oppenheimer approximation^ for the electronic 
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motion of H 2 + . 

The potential V, being the attraction of the two nuclei considered 
as fixed centers to the electron, is not spherically symmetric. 

As a result the total angular momentum of the electron is not con- 
served but its z-component is (the z-axis being defined as the line 
joining the two nuclei). Appropriate solutions of (3-6) can there- 
fore be written 


4 


L A <f 
e 




(3.7) 


where « and & are the polar angles of r (the vector coordinate of 
the electron). Substitution of (3*7) into (3 *6) yields 


J_ ^ \ + j t> — A VVJ 2 * £ E — V)1 ■= o j (3.8a) 

* iC-Ai** 2v - 


where /\ = 0, 1, 2, . . . correspond to the Y , IT J 

electronic states of H 2 + . The solutions can further be separated 
into even or odd electron parity ( i? tt — •& , ep ir + «p ) 

corresponding to gerade and ungerade classes. 


Our purpose is to derive '(3.8) from the radial equations of the 
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three-body system in the H limit. This provides not only a 
useful check of the radial equations themselves, but it allows for 
the verification of general rules which relate total angular momen- 
tum states of the three-body system to their Born-Oppenheimer counter- 
parts as one-body systems in an axially symmetric field. 

We shall do this in the following way: we shall actually show l* r 

the first couple of total jf radial equations in the H 2 + limit 
that they reduce to the form (3.8) for specific values of This 

will enable us to adduce a general relationship between total 

£ - states and the Born-Oppenheimer electronic states 
that they correspond to. Then we shall derive a formula relating 
the parity of the molecular electronic state to the parity and 
exchange character of the three -body wave function. 

A. / = 0 

From (2.2l) in the limit defined by (3-l) we have 


( E 

"B 


(3-9) 


where L* is given in (2.19). The latter is seen to depend 

differentially only on two coordinate: n and r 2 . This is already 

the seed of the Born-Oppenheimer limit which will be seen to apply 
to all l equations. 

Letting rj = |x ( + K + ) \ 


J 


f * 


o^wJL f — h. 



we find 


= i l - f. 1 1 + 2? + i L l + 2^ 

x n ^ J 1 3S »V *) *rp ? J»»- 3 S' 

(o *• 

Now letting r = Y| ,eos*> = ~ A I = *. We flnd further 

_ i ^ ^ + -J — 2_ £i’v> 2- 

*iv ~ >l v\ x tC~/S*fr ~iT> 'i* 

Thus (3.9) reduces to (3*8) for K =0. Ie. 0 

radial equation reduces to that for X states in the Born-Oppen- 
heimer limit. 

The difference between the Y. and higher A equations 
consists of terms proportional to (r 2 sin 2 tf ) 1 . We need only note 


the relation 
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so that the Born-Oppenhelmer equations in general can he written 


[C - 

r 


B. -C = 1 even parity 

For Z = 1 and *• = o, Eq. (2.2l) reduces to 


c e - v n ^ - ° 




4 *22.* C. E - lO J F,° — C> 


Comparing with Eq. (3.8b) we see that this is the Born-Oppenheimer 
equation for TT states. 


C. Z = 1 odd parity 

We get from Eq. (2.20) for ■£= 1, "H = 1 
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(3.8b) 
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(3.15a) 



It is also convenient to derive a redundant equation. Letting 


r x z* r 2 in ( 3 . 15 a) and using (2.18), we obtain 
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In order to show the Born-Oppenheimer limit of these equations 
we introduce the following transformation: 
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transformed of (3.15a) gives 

2 ^n + lw( £ E — l/) 

Subtracting gives 


(3.19b) 


Subtracting Eqs. (3.19 a ) and (3.19b), we get 

LU 1V + 2_v», e (_£■ - V) ] CH ' - Hj, ) 

* x 

which according to (3.8b) corresponds to H states in the Born- 
Oppenheimer limit ... 

Adding Eqs. (3.19a) and ( 3 . 19 b) we get 

[ C - <-4 
r 

which according to Eq. ( 3 . 8 b) corresponds 


*25 ( E - V) I («,' + Hi, ) = 


(3.21) 



r Li, - 2 + i«t(E-wi H H - 2 h . 1 = 

r * r 


H*. = o 
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(3.19a) 




22 


to n states in the Born -Oppenheimer limit. 

Therefore for ( = 1 we get two fl states and one E state. 

D. % = 2, even parity 

Eq (2.21) in this case involves the radial functions F° , F 2 , 
and F 2 . Making the transformation 

pO 1 ~ti 2 

h r= i ■. 

2 y 2 o 

c-2 • u 2 i a tJ 2 -ia 

F 2 = 1 H 2 e + H_ 2 e 

F 2 = H 2 e la + iH 2 e’ la 

2 2 “2 

yields three real equations for the functions H 2 , H 2 , H 2 . One can 

2 -2 0 

readily construct linear combinations of these equations which are 
the E, n, and A states of 

E. 1 = 2 , odd parity 

Here one makes the transformation 

F 1 =-iH 2 e la/2 + H 2 e ' lo/2 

2 i - 1 

~ = -H 2 e 1<s/2 + iH e‘ la/2 

F 1 1 “ 1 

2 

The resultant equations can then be shown in a similar manner to be 

+ 

equivalent to n and a states of . 

Therefore, it can be shown that for jf = 2 , there are E> Il> A 
even parity states and for odd parity there are n and A states. 


( 3 . 22 ) 

(3.23) 

(3.24) 


(3.25) 

( 3 . 26 ) 
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This checks the well-known rule that for an arbitrary { all 

£ 

A states except A = 0 are doubly degenerate (positive-negative) , 
and for £ states the parity of the l equation which yields this 
state is given by 

(-D 1 + K = i. 

If one considers the finite mass of the nuclei, then the degeneracy 
of the A f 0 states is lifted giving rise to the A " doubling 
phenomenon^. Our results also conform the rule^ {\ < f . 

In order to facilitate the discussion of the |=> - f*- -f» molecule 
we also give the relation of the electronic gerade and ungerade 
classes to the three-body symmetries. This relation is 

1E = 

which can readily be verified two ways. First both the left and right 
sides have precisely the same effect on the Hamiltonian, and second 
if one starts with a spatial configuration of the nuclei and the 
electron, then the operation of iE will yield the same final configu- 
ration as 

A specific example of the above is provided by ^ = 0 (there is 
only even parity in this case) and the l — 1 , odd parity cases. We 


(3.27) 


( 3 . 28 ) 


have shown above the radial equations in both cases approach the £ 
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equations of . Nevertheless the vector spherical harmonic 

portions of the complete wave functions are appropriate to 1 = 0 

and l - l respectively, and are hence different. The lowest 

X = o wave function corresponds to the ground £ state of H 

g 2 

whereas the lowest 1 = 1 odd parity wave function corresponds to 
the first excited rotational level of the electronic state of 
It is the latter state which in the ^ mesonic counterpart 

8 

of , namely the molecule p - y, - p, is primarily responsible 

for the absorption of ^ mesons in hydrogen. The original calcu- 

9 

lation which only corrects the Born-Oppenheimer approximation to 


the extent that the error is reduced from the order 


(m p./Mp)^^" to the order (m^/Mp) gives a decay rate of about 560 

-x 2 -i 

sec whereas the most recent experimental number is 464 + 42 sec . 

That the b ulk of the discrepancy is iri fact due to the Born-Oppenheimer 

approximation has been shown by a recent calculation of Wessel and 
3 

Fhillipson who attacked the problem as a three-body problem and 

obtained a value of 480 sec 1 . Their variational calculation did in 

fact omit certain coupling terms in the radial equation but these 

omissions are of the order of (m p. /Mp) 2 . (it should be noted that an 

accurate variational calculation has been made for this state in 

10 

which m 1 terms in the radial equations have been included . The 
calculation was not primarily intended for this application.) 
Furthermore, we see that the ground state, is gerade (iE = + l) 

and is derived from the t = 0 equation which is necessarily even 
parity ( * =0 ). From (3.23) it follows that £ 12 = + 1, so that 

the three-body spatial function is symmetric and its spin antisymmetri 
(singlet), i.e., a para state. On the other hand the capture takes 
place from an -l = 1 odd parity-state, the first rotational excited 
state of the T-g electronic state which according to (3.23) is 
S 12 - i.e. , a triplet or ortho state. Since the eigenvalue 
of (P is (-l) 


we can also write 
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APPENDIX A 


This appendix is concerned with the relation between the present 

Euler angles and another set of symmetric Euler angles that have 
11 11 

been introduced * . This particular alternate set differs only as 

concerns J . Specifically it is measured from the line of 
nodes to the larger of the moments of inertia in the particle plane. 
Letting ^ be the alternative azimuth, our objective is to find 
as a function 'f . We shall find that ^ depends 
additionally on r x , and © 12 , the significance of which we shall 
discuss further at the end of this appendix. 

The alternative definition of the azimuth is given in terms of 
a ratio of products of inertia in a coordinate system measured with 
respect to the center of mass. Letting fa be these coordinates of 
■the three particles , we find 


b - + b 



( A.l) 


fa * - j-( r. ** + ^ 
~ n 


(A.2) 


where > 1 = 1, 2, 3> are the masses of the three particles and 

p = i n • 

- \ 


(A.3) 
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(. ) - K i ) are the vectors from the third particle to each of 
the identical particles (although this analysis allows all three 
masses to he different. ) The components of the ^ are defined 
ty 


hr = L $i-x + ) 

where i' and j* are unit vectors along x'- and y 1 - axes in the 
particle plane. We have previously given, Eqs. (115) and (ll6), 

r a ^ ^ 

r]_ and r 2 in the i J ’ coordinate system. Substitution of these 
relations into (A.l) and (A.2) gives the following relations for ^ 

U*. - (l-hi ) (jf- j 01 *) _ |*„ K x J.W ii*- £ ) 

h »- 


Sly = -0 -hi) ’r, cosC4 - ±-6,*) * Kv 

^ 1 K 


Kk = - H' N - i6 (1 ) + ?ml4+ i 9 >') 


5*y - ill 'll Ctli t SP - - (]_-£»•} ’*•1 c&s( -£ 

u b 


(A.4) 


(A.5) 


(A.6) 


(A.7) 


(A.8) 
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Y 4 t* 


f J7 


Ht , cos i£'_iSti.) 




The definition of 


is given by 


13 


to.-* 1% - - 1 

1>X* ~ ^>YY 


where 
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Substituting (A.5) “ (A.IO) into (A.12) - (A.14) yields 
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Specializing to H 2 


limit, we have 
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where 
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C^x "V ) ®n. 


Using 


to*©,* 


^ *V *V|, — \t. 


we find 
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*\L **+*£) 


With the relation (A.16) one can then transform the wave 
12 

function 


c t >».k 

4 - t Hn(\, K t ( \. l4 ) ( e , § , ^ ) 

h«-A i 


into our form Ea. (2.1^), from which one derives the relation "between 
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(and hence ) and 

The radial equations which result are consequently precisely 
those that one would get directly from the formulae of Diehl et. al. ^ 
in the H 2 + limit. The y of their equations from the point 

of view of this limit stems from the fact that the axis to which 

is measured becomes identical with the line joining the nuclei, 
i.e., the z-axis, in the Bom-Oppenheimer approximation. 

In the case of all finite masses, however, there arises an 
additional group of terms in their equations. A comparison with our 
equations indicates that in general their equations are more compli- 
cated than ours. In particular some of the coefficients have a more 
complicated analytic behavior. We believe this is related to an ar- 
gument in I which stated that because our Euler angles depend only 
on the unit vectors r x and r^ and not on the magnitudes of r x and r 2 , 
we think our Euler angles are superior; for they retain the 
separation between angular and radial parts which was implicit in 
the expansion of the complete wave function itself (cf. (2.14) and 
(A.20)). The suggestion that we are making here is that a concrete 
manifestation of such a superiority may reside precisely in the 
differences in the analytic structure of the respective radial 
equations. 

The dependence of ^ on r^ and r 2 is explicitly exhibited in 

Eqs. (A. 16) and (A.17). In Appendix III of reference 1 we derived 

ih 2.5 

the relation between ours and the Hylleraas-Breit Euler angles ’ . 
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The latter have the disadvantage of not being symmetrical, but 

/■S A 

they too depend only on the unit vectors r^ and r 2 . Therefore 
the relations between the respective Euler angles are independent 
of r x and r 2 . According to the above argument we do not expect 
any significant difference in the analytic properties of Breit's 

*1 cr 

(P state) radial equations J as compared to our own. 
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Appendix B 

For three unequal mass particles, the kinetic energy in the center 
of mass system is 


T = 




(B.l) 


where 


f* , = w i K» 


( B . 2) 



+ M 


The Schrodinger equation in the c.m. 


system is thus 


[< + K + + ib (E-Wj { = ° • 

V (B.4) 

It has been shown how the individual' Laplacian may be written in terms 
of the Euler angles and residual coordinates X , 5 ; 9^* 
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where F x and F 2 are given in Eqs . (159, 162). Knowing then the effect 
of F x , F 2 and the cross term *7 K . on the vector spherical 

-v*. 

harmonics via their representation in terms of the A, , A^_ 

operators in ( 163, 64) and (2.13), one may readily derive the radial 
equations for the general case. As coupled equations in terms of the 
residual coordinates ,0 (2 , the radial equations are obtained 

by making the following substitution in ( 2 . 15 ): 
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Kv e,,. 2 0.o_ 


Since for m x m 2 the equations are no longer symmetric with 


(B.6) 


(B.7) 


(B.8) 


(B.9) 


respect to r x r 2 , it is necessary to solve the equations for all values of r x 
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and r 2 (as opposed to, say^ri y / r 2 and an appropriate boundary 
condition when mi is identical to m 2 ). By virtue of this there 
is no advantage in defining functions F^* and F* 
as they will also "be described by a coupled set of equations, 
rather than a single equation (2.2L) for the identic al particle 
case. For this reason ve give below the equations in terms of 
the coupled", fj functions, but involving the coordinates 
ri, r 2 , rxa. 
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where 


■*». - *- K* t 


' i , v£ 
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Mx 


Pt- -L X, . 
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'*' *■ **•* **<v 
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and 


^ | + 


(B.13) 


The resulting radial equations are invariant (except for a relative 

phase factor which does not affect the eigenvalue spectrum) under the 

simultaneous exchange of A, and ■**, ""V . However the equations 

are not (formal 1y) invariant under a simultaneous permutation of all 

three particles (i.e., simultaneous cyclic permutations of mi, life, 

m 3 and ri, r 2 , r^). This is apparently due to the fact that our 

choice of Euler angles singles out one particle, the particle with 

mass M, as the instantaneous origin. This in turn means that the 

line to which the azimuth f is measured depends asymmetrically on which 

particle we call the origin. As opposed to this the Euler angles of 
11 12 

Holmberg and Diehl et al. are invariant under the operation, and 
their radial equations are also. 

Since, however, the transformation between the respective radial 
functions can he worked out from Eq. (A.15), this additional symmetry 
can readily he recovered. It would seem that a more practical criterion 
of the utility of the equations is the analytic properties of the 
equations themselves. In this respect the discussion near the end of 
Appendix A may he more relevant. 
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TABLE I 

Coefficients of the Angular Derivatives in 
the Cross Term of the Kinetic Energy. a 
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aUie dependence on r x and r 2 is not included. All unincluded partial derivatives 
have zero coefficients. 











